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We present an ab initio analysis with density functional theory for superconductors (SCDFT) 
to understand the superconducting mechanism of doped layered nitrides /3-LirMNCl (M=Ti, Zr, 
and Hf). The current version of SCDFT is based on the Migdal-Eliashberg theory and has been 
shown to reproduce accurately experimental superconducting-transition temperatures T c of a wide 
range of phonon-mediated superconductors. In the present case, however, our calculated T c <4.3 K 
(M=Zr) and <10.5 K (M=Hf) are found to be less than a half of the experimental T c . In addition, 
T c obtained in the present calculation increases with the doping concentration x, opposite to that 
observed in the experiment. Our results indicate that we need to consider some elements missing in 
the present SCDFT based on the Migdal-Eliashberg theory. 

PACS numbers: 74.20.Pq, 74.25.Kc, 74.62.-c, 74.70.Dd 



I. INTRODUCTION 

Since the discovery of high-transition-tcmpcrature (T c ) 
superconductivity in cuprates,— superconductivity in lay- 
ered transition-metal compounds has been extensively 
studied. So far, a variety of exotic superconductors such 
as SrRu2C>4, Na^OoC^, iron arsenides and chalcogenides 
have been found. Among them, layered metallonitridc 
halides (MNX, M=Ti, Zr, Hf; X=C1, Br, I) discovered 
by Yamanaka et alW are of great interest, for which T c 
was found to be as high as 26 K at maximum! 3 ^ In fact, 
they had been the second highest record among transi- 
tion metal compounds until the recent discovery of iron- 
based superconductors. The mother compounds of MNX 
consist of stacking of two-dimensional MN block layers, 
sandwiched by intcrlaycr halogen ions. Since the inter- 
layer bonding is relatively weak, we can intercalate vari- 
ous Louis bases between the layers or deintercalate halo- 
gen atoms for the carrier doping. 4 Interestingly, while 
the mother compounds are non-magnetic band insula- 
tors, they become superconductors upon doping. This 
is in high contrast with many other superconducting lay- 
ered transition metal compounds, for which various mag- 
netic phases reside in the vicinity of the superconducting 
phase. 

When we consider the mechanism of superconductiv- 
ity in the doped MNX, we have two possibilities. One is 
a conventional scenario based on the Migdal-Eliashberg 
(ME) theorypHS where the superconductivity is medi- 
ated by phonons, and the damping and retardation ef- 
fects are treated by self-consistent perturbation theory, 
retaining the lowest-order dressed-phonon and dressed- 
Coulomb contributions to the self energy. The momen- 
tum dependences of the gap function, self energy, and 
screened Coulomb interaction is neglected within the en- 
ergy scale of phonons. In the ME theory, the screened 
Coulomb interaction is conventionally addressed with the 
random-phase approximation (RPA) and the frequency 



dependence is totally ignored. The other scenario is to 
consider the factors neglected in the ME theory, in which 
we consider various pairing glues other than phonons, 
such as spin, charge, or orbital fluctuations. In fact, 
for some layered transition metal superconductors, it has 
been believed that the pairing mechanism is unconven- 
tional. For example, the anisotropic d-wave pairing in 
cuprates and triplet pairing in Sr 2 Ru04 arc difficult to 
understand within the ME theory. As in the case of iron- 
based superconductors, when the electron-phonon cou- 
pling is weakp we can also safely exclude the possibility 
of the conventional scenario. On the other hand, the 
situation is not so simple for doped MNX. 

Experimentally, the Knight shift decreases toward 
zero b elow T c , indicating that the Cooper pair is spin 
singlet E2LLU The break-junction and scanning tunneling 
spectroscopy measurements sug gest that the gap func- 
tion does not have nodesP^HSl While these are typical 
behaviors of conventional superconductors, there are also 
several observations which support unconventional sce- 
nario. Measurements of uniform spin susceptibility^ 
and specific heaps indicate that the electronic density 
of states (DOS) at the Fermi energy (Ep) is small and 
the estimated total electron phonon coupling A is smaller 
than 0.22. On the other hand, the specific-healP^ and 
tunncling-current measurementaI2HHl give a rather large 
superconducting gap ratio ^^~4.6-5.2 or even larger, 
implying the strong effective pairing interaction between 
electrons. Also, the nitrogen isotope effect coefficient has 
been reported to be markedly small as ~0.07p2E2 com _ 
pared with the value of 0.5 for the BCS superconductor.^ 
Doping dependence of T c (Ref. ITM2TJ) is also difficult to 
understand in terms of the conventional scenario; par- 
ticularly in /3-Lia.ZrNCl, it becomes higher with decreas- 
ing doping and exhibits the highest value at the bor- 
der of the Anderson insulating phase and the super- 
conducting phase. On top of these, more recently, it 
has been reported that the spin-lattice relaxation rate 



of the nuclear- magne tic resonance does not exhibit the 
coherence peakpEU anc [ muon spin rotation measure- 
ment indicates that the gap function has a significant 
anisotropyP2 

As for theoretical studies, both the phonon and elec- 
tron mechanisms have been examined and proposed. For 
the electron mechanism, a possibility of spin-fluctuation- 
mediated superconduct ivity has been studied within the 
lattice Fermion modelpiE2H2U For the phonon mecha- 
nism, ab initio analyses based on the ME theory have 
widely been performed: Heid and BohnerP^ calculated 
the phonon spectrum and the electron-phonon coupling 
using density functional perturbation theory and esti- 
mated T c using the dirty-limit gap equationPB They 
found that A is as large as 0.5 for /3-LijygZrNCl and 
the value is too small to explain the experimental value, 
T c — 16 K, with setting the Coulomb pseudo potential fj,* 
(Ref.inU) to a typical value of -0.10. Weht et alW^ evalu- 
ated the Hopficld-McMillan parameter 77, also a measure 
of the electron lattice coupling, for /3-Na .25HfNCl using 
the Gaspari-Gyorffy formula.^ They found that rynf— 0.5 
eV/A 2 and t?n— 0.2 eV/A 2 , which are one order of mag- 
nitude smaller than those of the binary transition-metal- 
nitride superconductors with T c ~20 KP2 On the other 
hand, recently, Yin et al. 35 proposed that A is estimated 
to be larger if we employ the hybrid Heyd-Scuseria- 
Ernzerhof functional and T c giv en by the Macmillan- 
Allen-Dines (MAD) formulaPE^ with ^4* =0.1 shows a 
nice agreement with experiments. It should be noted 
here that there are many examples where fj,* substan- 
tially (larger than 50%) deviates from the typical value 
~ O.lpSso the problem whether the experimental T c can 
be described within the ME theory is yet to be clear. We 
also note that there have been several studies which try 
to consider the physics beyond the ME theory; for exam- 
ple, Bill et al. showed that acoustic plasmon can enhance 
the phonon- mediated superconductivityPB There are also 
several studies which propose that disorder can enhance 
the pairing instability in sy stems with phonon-mediated 
attractive interactions! 41 * 42 ! 

For the electronic interaction, there are many efforts 
for describing it from first principles. Lee et alW^ pro- 
posed a method to evaluate the bare Coulomb parame- 
ter fi with ab initio RPA calculations. However, there 
still remains an empirical parameter representing the 
electronic energy scale, needed for obtaining the renor- 
malized value /i* . On the other hand, Liiders et al. 
recently developed a parameter-free method to calcu- 
late T c based on de nsity functional theory for super- 
conductors (SCDFT) j 44 * 45 ! where electron-phonon cou- 
pling and electron-electron interactions are nonempiri- 
cally treated. Since the exchange-correlation kernels are 
constructed on the basis of the ME theory, the current- 
stage SCDFT has a close correspondence with the ME 
theory. This method has been applied to a wide range of 
typical phonon-mediated superconductors such as simple 
metals, 46 MgB 2 j !47 calcium-intercalated graphitep^ Li, K 
and Al in high pressure®! and CaBeSip^ and the relia- 



bility on the quantitative aspect has been examined and 
established, where T c s have shown agreements with ex- 
perimental results within a range of few K. 

In this work, we performed an SCDFT analysis for the 
lithium-doped /3-MNC1 (M=Ti, Zr, Hf) superconductors 
to examine whether they can be described by the cur- 
rent SCDFT formalism based on the ME theory. We 
calculate electronic structures, phonon spectra, electron- 
phonon couplings, screened Coulomb interactions in the 
RPA level to estimate superconducting T c and examine 
the results focusing on the metallic-atom and doping de- 
pendences. We found that T c estimated by SCDFT is 
at maximum half of the experimental T c and its doping 
dependence is opposite to experiments, both of which im- 
ply that, as in the unconventional superconductors, the 
present superconductors require some elements missing 
in the conventional ME theory. 

In Sec. [nj we review the SCDFT formalisrri 4 ^ 6 ^ and 
describe our computational detail. Applications to alu- 
minum and niobium are performed in Sec. |III| to test the 
reliability of our calculations. We show results for /3- 
Lij-MNCl in Sec. |IV| Section[V]is devoted to discussions. 
Summary and outlook are drawn in Sec. |VI| 



II. METHOD 



In SCDFT, we solve the following gap equatio n 1 44 * 46 ! 



^nk — ^-"rikh^iTk. n / J *^iikn'k.' 



nV 



tan h[(3/2)E», 



] A*.(1) 



Here, n and k denote the band index and crystal mo- 
mentum, respectively, A is the gap function, and /3 is 
the inverse temperature. The energy E n \^ is defined as 
E n k= vHk + ^nk an d £«k = e «k - ^ is tne one-electron 
energy measured from the chemical potential /1, where 
e„k is obtained by solving the normal Kohn-Sham equa- 
tion in density functional theory (DFT) 



"HKsl^nk) = Enkl^nk) 



(2) 



with "Hks an d \<Pnk) being the Kohn-Sham Hamiltonian 
and the Bloch state, respectively. In Eq. ([!]), Z and K are 
the exchange-correlation kernels in the SCDFT formal- 
ism. For the superconducting phase, the self-consistent 
solution of Eq.([l]) becomes nonzero; thereby we deter- 
mine T c . 

The exchange-correlation kernels Z and K, are de- 
fined by the second-order derivative of the exchange- 
correlation energy with respect to the electronic anoma- 
lous density. In the present study, we consider the 
three terms^ shown in Fig. u\ to derive Z ph [(a)], K. ph 
[(b)], and K. cl [(c)l Here, Z ph (/C ph ) is the lowest- 
order contribution^^ to the exchange-correlation ker- 
nel due to the interaction between phonons and normal 
(anomalous) electrons, and /C el represents the exchange- 
correlation interaction between anomalous electrons via 
screened Coulomb interaction in the low- frequency limit. 
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FIG. 1: Schematic pictures of the exchange-correlation en- 
ergy to derive the three kernels, Z ph (a), fC ph (b), and KL cl 
(c). The solid line with arrows in the same (opposite) direc- 
tion denotes electronic normal (anomalous) propagator. The 
dashed line denotes phononic propagator, and the bold wavy 
line denotes the screened electronic Coulomb interaction in 
the low- frequency limit, respectively. 



These kernels are given in the form averaged over the 
phonon energy ui and the electronic energy £' ranging 
from — fi to oo as 



Z^Z^ nk ) 



tanh[G0/2)& 



r 

dg I dua 2 F{uj) 



x[J(t*k,?,u) + J($,*,-€,w)] 



(3) 



and 



^nk,n'k' _ ^ (snkiCn'k') 



tanh[(/3/2)4 k ]tanh[(/3/2)^, k ,]iV(0) 
x / dwa 2 F(u})[l(Uk,C>k>,u) 

-'&*,-&*,«)], (4) 

respectively, where N(0) is the electronic density of states 
(DOS) at the Fermi energy per spin and a 2 F(ui) is the 
Eliashberg function defined as 

Here, w q „ is the phonon frequency of the mode q in the 
branch v. Hereafter, we use k and q to specify the mo- 
mentum of electrons and phonons, respectively. The 7 q „ 
is the phonon linewidth and given as 

7 q „ = 2irUq V Yl l5k+qn' 1 kn| 2 ^(&+qn')<*(£kn) (6) 

knn' 

with 3k+qn' kn being the electron-lattice coupling as 

5k+qn',kn = (<Pn'k+q|<% I /|Pnk)- (7) 

Here, SV^y is the potential derivative with respect to qy 
phonon normal mode, which couples the two electronic 
states \ipnk) and |<p n 'k+q)- 



The integrand function I in Eq. Q is written as 



,fc 



*?(£'+<•>) 



oK' 



/3(?-K>) 



i-a-u 



c-e'+w 



(8) 



where fp and np are the Fermionic and Bosonic distri- 
bution functions, respectively. The integrand function J 
in Eq. Q is given as 



■7&£,w)=jte£>)-jter,-«) 



(9) 



with 



J(Z,?,u)-- 



MQ + npjw) 
t-?-w 



mo-mc-oj) 



e - e - « 



-0MZ-u)M-S + o>) 



(10) 



The kernel /C o1 ascribed to the screened Coulomb in- 
teraction between the Cooper pairs is given as 



K. 



„k,/k-- / dr j dr'cp* nk (r)(p* n _ k (r')W{r,r') 

XVn'k'fT^n'-k'fr'), (11) 



where W(r,r') is the static screened Coulomb interac- 
tion. In this paper we evaluate this term with two kinds 
of approximations, i.e., the Thomas- Fermi approxima- 
tion (TFA) (Ref. HS] ) and the random-phase approxima- 
tion (RPA)PEUHI Within TFA, the screened Coulomb 
interaction is written with the Thomas-Fermi wavenum- 
ber &tfa — \/8TriV(0) as 



WTFA(r,r') = 



-fe T PA|r— p'| 



(12) 



Since Wtfa(i", r') is given by the sum of its Fourier com- 
ponents as 



VF T FA(r,r')=5Z 



-in 



k,G 



k+GI 2 + fc2 



i(k+G)-(r-r': 



,(13) 



TFA 



substitution of Eq. (1121) into Eq.(ll) yields 

1 



jj-el.TFA _ 47T ^-s 



n'k' 



<) ^ |k-k'+G| 2 +fc| FA 



Pl£ (G)f (14) 



with 



n'k' 
Pnk 



;G)= f dr^^CrJ^^-^^-VnkW, (15) 
Jn 



where G is the reciprocal lattice vector, fi is the volume 
of the unit cell, and J Q denotes the integration within the 
unit cell. In this approximation, the screened interaction 
is isotropic and the screening length is determined by the 
electronic DOS at the Fermi energy. This interaction is 
nothing but that derived from the RPA in the vacuum 



and in the long- wavelength limit, where the local field 
effect of crystal is totally ignored in the calculation of 
screening. 

We next consider the RPA expression of the screened 
Coulomb interaction in the low-frequency limit. In gen- 
eral, the screened Coulomb interaction of crystal is ex- 
pressed with the double-Fourier transform as 54 

W(r,r')=J2 e i ( K+G )- r ^ GG '(K) e - i ( K + G '> r ' (16) 



KGG' 



with 



IF GG ,(K)=4^- 



K + GI 



£ GG 



,(K) 



K + G' 



(17) 



Here, K is the wavevector in the first Brillouin zone 
and £gg'(K) is the symmetrized dielectric matrixpS 
which is calculated using the following expression 



£gg'(K)=6 G g' -47T- 



1 



K + G 



■Xgg'(K) 



1 



K + G 



7T(18) 



with x° being the independent-particle polarization as 



XGG'(K)-- 2^ " 



n,n ; ,k 



£nk — £n'k+K 



[{p: k k+K (G)}>l k+K (G') 



(19) 



After inserting Eq. \1Q\ into Eq. ( 11 ) and with K=k k', 
we obtain 



/c: 



, RPA _ 4tt ^pl%, (G)e GG , (k-kO{^ k , (GO}* 



GG' 



|k-k' + G||k-k' 



(20) 



The resulting kernel /C el ' RPA includes the local-field effect 
on the screening in the real crystal. In this expression, 
£el,HPA has a s i n g U i ar i ty at k=k' and G=0 or G'=0. 
This singularity is avoided following Ref . 1551 

Figure [2] describes the flowchart of the present SCDFT 
calculations. Starting from solving the Kolm-Sham equa- 
tion in Eq. (21), we next perform density-functional per- 
turbation calculations to obtain a 2 F(ui), which generates 
electron-phonon kernels Z? h [Eq. ^} and /C ph [Eq. Q]. 
Also, we calculate the electron-electron kernel /C cl for 
TFA with Eq. ([14} and for RPA with Eq. §20§. With 
these information as inputs, the gap equation in Eq. M 
is solved. According to the two types of /C el ' TFA and 
/C el ' RPA , two values of T c are estimated, which are de- 
noted as ySCDFT-TFA and ySCDFT-RPA p rom compari - 

son between T C SCDFT - TFA and r c SCDFT " RPA , we study the 
effects beyond the uniform and local electronic screening 
on the transition temperature. 

We remark some technical details when we solve the 
gap equation in Eq.Q. In the present calculation, the 
states in the gap equation, labeled by nk, are generated 
by sampling; they are generated so that their density 
distributes logarithmically in the energy range of inter- 
est. For this purpose, with a given band n, the Brillouin 
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FIG. 2: Procedure to obtain the superconducting transition 
temperatures through density functional theory for supercon- 
ductors. 



zone (BZ) is divided into four regions, following a cri- 
terion for k; (i) |f nk | < uj , (ii) cj < |£„k| < Wi, (iii) 
wi < |£„k| < w 2 , (iv) u> 2 < |£„k|, where £ nk is the energy 
of the sampling state and uii {i—Q, 1, 2) are the energy cri- 
teria. In the present case, ujq, u\, and uj-i are set to 0.002 
eV, 0.02 eV, and 0.35 eV, respectively. Next, to realize 
logarithmic distribution of the sampled states, we intro- 
duce an acceptance ratio p for each region; p=\ for (i), 
p=0.1 for (ii), p=0.01 for (iii) and p=0.002 for (iv). With 
this acceptance ratio, the randomly-generated sampling 
states form a logarithmic distribution. For the resulting 
nk points, the energy £„k and the exchange-correlation 
kernels Z n \^ and /C ra k,n'k' are ev aluat ed with the linear 
tetrahedron interpolation schemeS^EU for the original ab 
initio data. On the ground that the density of the sam- 
pling points is not uniform, the gap equation in Eq.(l) is 
rewritten as 



A, 






>/Cr, 



tanh[(/3/2)£;. n 

-Era'k' 



(21) 



with W n k=V na /N na being the weight normalized as 
Sk W / nk=l for each n. Here, a specifies the region where 
the sampling state nk belongs, and V na is the volume of 
the region a in the first BZ for band n, estimated by 
the frequency of sampling.^ The total number of the ac- 
cepted sampling points, N na , also depends on a and n. 
We found that 10,000 sampling states per band cross- 
ing the Fermi energy assure a convergence within a few 
percent. 

As a reference, we also calculate T c using the 
McMillan- Allen-Dynes formula (MAD) 



T c , 



Win. 

1.2 



exp 



1.04(1 + A) 



A- M *(1 + 0.62A) 



(22) 



The parameters A and wi n are defined using the Eliash- 
berg function a 2 F as 



A = 2 / duj 



a 2 F{uj) 



and 



wi n = exp 



jdu 



U! 



a 2 F(uj) 



lncj 



Jdu 



a 2 F(u) 



(23) 



(24) 



respectively. The parameter fi* in Eq. ( 22 ) is an effective 



Coulomb pseudopotential in the narrow energy region 
within the Debye frequency wd above/below the Fermi 
level. Within RPA, 43 /i* is written with the renormaliza- 
tion formulaPD as 



j« = 



/'- 



1 + iAu[ 



(25) 



where E c \ is a parameter defining the electronic energy 
scale. It is conventionally set to the Fermi energy mea- 
sured from the conduction-band bottom, assuming that 
the band dispersion is parabolic. The bare fj, is defined 

as 



,i = A(0)((/C el » FS , 



(26) 



where ((/C o1 ))fs is the average of the electronic interaction 
over the Fermi surface and is defined as 



((K el »FS 



1 



A(0)< 



E *(&.kW&.'kO£k 1 n'k'- (27) 



Since [i is estimated for both of /C cl < TFA [Eq. (|14|] 
and /C el,RPA [Eq. (20)], we distinguish them as /itfa 



and /xrpa- Thus, we derive the two T c values, i.e., 
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FIG. 3: Procedure to obtain the superconducting transition 
temperatures based on the McMillan-Allen-Dynes formula. 



T c MAE Vtfa) and ^c MAD (/4pa)- We summarize in FigJ3J 
the procedure of the MAD calculations. 

In the MAD analysis, the quantitative reliability of the 
resulting T c is somewhat ambiguous, because the renor- 
malization formula in Eq. ( 25 1 includes an empirical pa- 
rameter E c \. An advantage of the SCDFT framework 
over the MAD analysis is that the former does not contain 
any empirical parameters such as E c \. In the present pa- 
per, we study quantitatively the electronic-interaction ef- 
fect on T c with systematic comparisons between SCDFT 
and MAD, and between TFA and RPA. 



III. SIMPLE METALS 

To check the reliability of our SCDFT code, we ap- 
plied it to benchmark systems; aluminum and niobium, 
which are weak-coupling and strong-coupling BCS super- 
conductors, respectively. The comparison between our 
results and experiments, together with other SCDFT re- 
sult based on TFA,E£lare shown in Table[l] The table also 
contains the T c values based on the MAD formula and 
several key parameters used in the MAD calculations. 

Our calculated T C SCDFT - TFA and T C SCDFT " RPA show 
a good agreement with the previous calculation 
( T c S ?c? FT " TFA ) and the experiments (T^). The small 
deviation between the present and previous calculations 
comes from minor differences in a 2 F used for the calcu- 
lation of the phononic kernels. It is interesting to note 



TABLE I: Our calculated superconducting transition tem- 
perature T c based on density functional theory for super- 
conductors (SCDFT), together with the result bas ed on the 
McMillan- Allen-Dynes (MAD) formula in Eq.(22]). Exper- 
imental and previous SCDFT results are listed for compari- 
son. Parameters used in the MAD calculations are also given. 
The jUrpa and /irpa parameters are compared with the other 
theoretical results. 
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"SCDFT result with /C cl 
from RofHol 

Experimental values taken from Ref l61l 
c Ab initio RPA result in RefB3l 
d Ab initio RPA result in Refl60l 
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that the MAD formula significantly overestimates T c [see 
T c MAD (a4fa) and T c MAE Vrpa) in the table]. This error 
can partly be attributed to the oversimplification of the 
elect ron ic-interaction kernel in MAD, i.e., the use of \x in 
Eq. (26), 62 ignoring the wavenumber dependence of the 
screened Coulomb interaction. In this approximation, 
the Coulomb interaction is approximated by a 5 function 
in the real space, so that electrons forming the Cooper 
pair do not feel distant Coulomb repulsion. As a result, 
the electronic interaction is estimated to be smaller in 
MAD than SCDFT considering the wavenumber depen- 
dence of the Coulomb interaction and T c obtained by 
the MAD formula tends to be higher than T c estimated 
within SCDFT. We also note that the difference between 
TFA and RPA is small, which indicates that the local- 
field effect is not so significant for Al and Nb. 

We note that the deviation of our j , c SCDFT from the 
experiment for Al is relatively large if we consider the 
T c ratio. As suggested in Ref 46, in the superconductors 
with small gap with T c below a few K, the T c value is 
a consequence of a subtle balance between the electron- 
phonon coupling and the screened electron-electron in- 
teraction, thereby systematic errors depending on the 
computational details of approximations are strongly en- 
hanced. This type of error arising from a subtraction of 
two nearly equal quantities does not occur in the super- 
conductors with relatively high T c and do not need to be 
worried in the following sections. 



IV. LAYERED NITRIDE CHLORIDES 
A. Computational Detail 

Ab initio electronic and lattice-dynamical calculations 
were performed for lithium-doped layered metallonitride 
chrolides, /3-Li^MNCl (M=Ti, Zr, and Hf) with the 
three doping rate x=0.0, 0.3, and 0.5 with Quantum 
Espresso package,— where the local density approxi- 
mation (LDA) with the parameterization by Perdew and 
ZungerS3E3 and the Troullier-Martins norm- conserving 
pscudopotentiala^were employed. The pscudopotentials 
for Ti, Zr, and Hf were generated in the semicore configu- 
rations of (3s) 20 (3p) 60 (3d) 20 , (4s) 20 (4p) 60 (4d) 20 , and 
(5s) 20 (5p) 60 (5(i) 20 , respectively. The scalar-relativistic 
correction^ was applied to the Zr and Hf pseudopoten- 
tials. The Li pseudopotentials were supplemented with 
the partial core correction.^ The lattice parameters and 
internal coordinates were fully optimized. We checked 
that our calculations for the relaxed structures and the 
T-point phonon frequencies well reproduce the experi- 
ments (Appendix |A[) . 

The cutoff energies in the wavefunctions were set to 
75 Ry, 80 Ry, and 90 Ry in the Ti, Zr, and Hf com- 
pounds, respectively. The charge density was described 
with the 8x8x8 k points in the Monkhorst-Pack grids 69 
and the electronic DOS was calculated using the denser 
32x32x32 k grids. Phonon dynamical matrices were 



calculated with the 4x4x4 q points, while the electron- 
phonon couplings {<?k+qn' kni m Eq- y\ were calculated 
with the 4x4x4 q phonon modes and the 32x32x32 k 
wavefunctions. The polarization matrix (xq G ;(K)} in 
Eq. ( 19 ) was expanded with the plane waves with the 



cutoff of 12.8 Ry and calculated on the 5x5x5 k grid. 
The Kohn-Sham states within a range [—12 eV, 36 eV] 
were included in the polarization calculation. The Fermi- 
surface integrals involving 7 q „ in Eq. |6|) and fj, in Eq. ( 26 ) 



were performed using the first-order and zeroth-order 
of the Hermite-Gaussian smearingpS respectively. The 
smearing width was set to 0.020 Ry in x=0.5 and 0.025 
Ry in £=0.3. The SCDFT calculation in Eq. @ was 
done with 10,000 sampling k points for bands crossing 
the Fermi energy and 500 points for the others, where we 
considered 25 valence and 25 conduction bands. 



B. Atomic geometry 

We show in Fig. El atomic structures of /3-ZrNCl (a) 
and /3-Li .5ZrNCl (b). We display the side view along 
the a axis. The atomic structure consists of stacking 
block layers made from Zr (gray), N (blue), and CI (red) 
atoms. On inserting Li (yellow) in between the block 
layers, the crystal slightly expands. Here the positions 
of the Zr, N, and CI atoms, as well as the intercalated 
Li atom have been determined with experiment s.EBESl j^ 
is also experimentally confirmed that the stacking order 
of the layers for doped systems is different from that of 
the mother compounds P^ Th e Hf compounds basically 
take the same structure! 76 * 77 ^ The details of the atomic 
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FIG. 4: (Color online) Atomic geometry of /3-ZrNCl (a) and /3- 
Li^ZrNCl (b) along the a axis, displayed with the conventional 
cell (solid line). Gray, blue, red, and yellow spheres stand for 
Zr, N, CI, and Li atoms, respectively. The space group is R3m 
(No. 166). The Li atom is located at the 3a site in terms of 
the Wyckoff position, whereas the other atoms are located at 
the 6c sites, whic h are specified by one parameter z (see also 
Table. P and [!!l]|. 



TABLE II: Our calculated structural parameters for undoped 
/3-MNC1. 



M 


Ti 




Zr 




Hf 




present 


present 


expt. a 


theory" 


present expt. a 


afAl 


3.342 


3.553 


3.605 


3.556 


3.488 3.577 


cfAl 


26.164 


26.879 


27.672 


27.178 


26.762 27.711 


z(M) 


0.120 


0.119 


0.119 


0.118 


0.119 0.120 


z(N) 


0.197 


0.198 


0.198 


0.198 


0.198 0.198 


z(CV) 


0.388 


0.386 


0.388 


0.386 


0.387 0.388 



a Single-crystal X-ray diffraction measurement 

taken from Ref. 1741 

b Ab initio DFT-LDA full-potential result in Ref . 1751 



TABLE III: Our calculated structural parameters for doped 
/3-LL.MNCL 



M 


Ti 


Zr 




Hf 






present 


present 


cxpt. a 


present 


expt." 


X 


0.3 0.5 


0.3 0.5 


0.20 


0.3 0.5 


0.29 


a\A] 


3.394 3.406 


3.604 3.616 


3.591 


3.541 3.561 


3.589 


cfAl 


27.740 27.524 


28.163 27.942 


27.839 


28.173 27.795 


29.722 


z(M) 


0.207 0.206 


0.209 0.208 


0.213 


0.208 0.207 


0.208 


a(N) 


0.136 0.134 


0.135 0.132 


0.136 


0.135 0.133 


0.137 


2 (C1) 


0.391 0.385 


0.388 0.383 


0.389 


0.389 0.384 


0.394 



"Single-crystal X-ray diffraction measurement in Ref. 1731 
Powder neutron diffraction measurement for /3-Na^HfNCl in Ref. 1771 



C. Electronic Bandstructure 

We show in Fig. [5] our calculated band structure and 
electronic DOS of /3-Li^MNCl (x=0, 0.3, and 0.5; M=Ti, 
Zr, and Hf) for the range of [—5 eV, 4 eV]. The undoped 
system {x—0) is an insulator with a band gap of ^0.5 
eV for M=Ti and -2 eV for M=Zr and Hf. Upon doping 
(from the left to right panels) , we see that carriers are ac- 
commodated in the conduction bottom, without chang- 
ing its two-dimensional dispersion significantly. While 
the DOS of this two-dimensional band is almost constant 



as a function of energy, N(0) increases (see Table IV I by 
carrier doping. For the material dependence (from the 
top to bottom panels), we see that the band structures 
have several common features; for example, the conduc- 
tion bottom always locates at the K point. On the other 
hand, DOS around the Fermi level is larger for the lighter 
element systems. In particular, the Ti compounds have 
high N(0) because the van-Hove singularity is closer to 
the conduction bottom, which seems favorable to have a 
strong electron-phonon coupling. 



D. Phonon Spectra 



configuration are found in Ref. |U 

We summarize in Table III] the optimized structural pa- 
rameters for the non-doped system /3-MNC1. The lattice 
parameters are approximately 3% smaller than the ex- 
periment and the internal z parameters agree well with 
it. Although the Ti compound with the /3 structure is 
not synthesized experimentally, in the present paper, we 
study this material as a reference and compare it with 
other materials. The optimized internal parameters of 
the Ti compound are very similar to those of the Zr and 
Hf compounds, whereas the lattice parameters are some- 
what smaller, which is reasonably understood in terms of 
smallncss of the ionic radius of Ti. 



We next show in Table IIIII our calculated structural 
parameters for the doped system /S-Li^MNCl. For the Zr 
compounds, the theoretical parameters agree well with 
experiments. For the Hf compounds, the experimental c 
parameter is somewhat larger than the theoretical value, 
which may originate from the difference in the interca- 
lated atoms, Li and Na. For the reliable description of 
the atomic structure of the incommensurate filling x=0.3, 
the virtual crystal approximation that the Li ionic charge 
is artificially reduced from the nominal value of 3 to 2.6 
was employed. This is critical to reproduce the experi- 
ment; a conventional treatment of electron doping with a 
uniform compensating positive background charge leads 
to an inaccurate description for the lattice parameters, 
a deviation from the experiment by ~10 %. In our vir- 
tual crystal approximation, the error is ~1% for the Zr 
compound. 



We next show the phonon dispersions in Fig. [6] The 
doping (material) dependence can be seen by compar- 
ing the panels from left to right (from top to bottom). 
The total spectrum consists of two bunches of bands; 
the high-frequency bands around 500-700 cm -1 and the 
low- frequency bands around 0-400 cm . The high- 
frequency bands mainly originate from the vibrations of 
nitrogen atom, while the low-frequency ones are formed 



TABLE IV: Density of states N(0) at the Fermi energy [/(eV 
spin f.u.)] a nd parameters in the McMillan-Allen-Dynes for- 
mula in Eq.(22i. For V(0), A, and wi n , our calculated results 



are compared with the experimental and previous theoretical 
results. 



M 



Ti 



Zr 



Hf 



X 


0.3 0.5 


0.3 0.5 


0.3 0.5 


N(0) 


0.342 1.675 


0.191 0.662 


0.176 0.437 


iV«pt.(0) 


- 


0.23° 


0.25 6 


AT LDA (0) 


- 


0.169-0.146 c 


0.19 d 


A 


0.286 0.411 


0.552 0.982 


0.827 1.292 


Acxpt. 


- 


0.22" 


- 


Alda 


- 


0.521-0.508 c 


- 


Wib [K] 


459 429 


421 305 


324 204 


wlda [K] 


- 


422-415 c 


- 


fJ-TFA 


0.771 1.705 


0.329 0.576 


0.413 0.367 


MRPA 


0.238 1.066 


0.161 0.441 


0.161 0.263 


Mtfa 


0.281 0.331 


0.177 0.218 


0.200 0.178 


Mrpa 


0.155 0.296 


0.113 0.196 


0.114 0.149 



Spccific-hcat measurement taken from Ref |16l 
Magnetic susceptibility measurement taken from RefjJ5] 
c Ab initio DFT-LDA result for ^-Lii/gZrNCl in Reff29l 
d Ab initio DFT-LDA result for ,3-Nao. 2 oHfNCl in Ref[32] 
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FIG. 5: (Color online) Our calculated band structures and density of states (DOS) of /3-LizMNCl (x=0, 0.3, and 0.5; M=Ti, 
Zr, and Hf). The DOS is given in the unit of eV~ 1 spin _1 (f.u.) _1 . The DOS was calculated with the tetrahedron integration 
method^ for the electronic bands interpolated onto 32x32x32 k mesh. Dashed lines represent the Fermi energy. We employed 
a conventional set of special points following Ref. [29j The points Kz and Mz stand for the K+Z and M+Z points, respectively. 



by the vibrations of lithium, chlorine, and transition- 
metal atoms. For the high-frequency bands, Li doping 
causes appreciable band narrowing and a red shift. On 
the other hand, for the low-frequency bands, carrier dop- 
ing basically does not change the overall profile so drasti- 
cally, but some characteristic change is observed for the 
Hf compounds; along the Z-Kz line, significant mode 
softening near 0-100 cm -1 occurs. 

Figure [7] displays the phonon DOS, where the upper 
three panels describe the doping dependence in the same 
material and the lower three panels compare the results 
of the three materials with the same doping rate. The 
figures show that the carrier doping and the substitu- 
tion of light elements with heavy elements lead to a red 
shift of the phonon DOS, especially for the low-frequency 
range of 0-200 cm" 1 . The origin is explained as fol- 
lows: For the effect of the carrier doping, N(0) becomes 
larger (Table IV I, so that the electronic polarization is 
enhanced. Thus electronic screening works more effec- 
tively and weakens the atomic interaction. On the other 
hand, heavier atoms vibrate more slowly, so that they 



generally have lower phonon frequencies. 



E. Electron-Phonon Coupling 



To estimate the strength of the electron-phonon cou- 



pling, we calculated the mode-dependent A qi , as 



1 



Iqu 



V " 7dV(0) w, 2 



(28) 



q// 



Note that A q „ is finite for metals because 7 q „ in the 
numerator is defined via the Fermi-surface integral [see 
Eq. J6b], The results are displayed in Fig. [8] We also 
plot in the right side the Eliashberg function a 2 F(uj) in 
Eq. ([5]) and accumulated-frequency function X(uj) defined 
as 



\(u) = 2 / du 

Jo 



,a 2 F{uj') 



W 



(29) 



The strength of the electron-phonon coupling can be 
measured by A(w) and the MAD parameter A in Eq. (23) 
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FIG. 6: (Color online) Phonon dispersions of /3-Li^MNCl (x=0, 0.3, and 0.5; M=Ti, Zr, and Hf). The acoustic sum rukP'was 
applied for the T-point modes. 
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FIG. 7: (Color online) (a)-(c): Doping dependence of phonon density of states of /3-LizMNCl for each transition metal, (d)-(f): 
Transition-metal dependence of phonon density of states in /?-Li x MNCl with a fixed doping ratio. The density of states are 
calculated with the tetrahedron integration method*^ for the phonon branches interpolated with a 32x32x32 fine mesh. 



is defined as the uj— >oo limit of A(w). We see that carrier doping enhances the electron-phonon coupling (compare 
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FIG. 8: (Color online) Electron-phonon coupli ngs for /3-Li^MNCl (x=0.3 and 0.5; M=Ti, Zr, and Hf); the mode-dependent 
electron-phonon coupling coefficient A q „ in Eq. (28 1 (left pane l) , Eliashberg spectral function a 2 F(u>) in Eq. (jsl) (right panel), 
and frequency-dependent coupling coefficient A(cu) in Eq. (29 1 [(blue) solid curves in the right panel]. The Eliashberg function 
is calculated with the tetrahedron integration method^ for the phonon branches and linewidths interpolated with a 32x32x32 
fine mesh. 



the left and right panels). Similarly, the electron-phonon 
coupling is strong in the systems with heavier transition- 
metal elements (compare the top, middle, and bottom 
panels) . 



Table |IV| lists our calculated A in Eq. (231 to quantify 



the electron-phonon coupling of each material and each 
doping rate. In terms of A, we expect high T c for heavier 
transition-metal compounds in the high-doping regime. 
However, on the other hand, the averaged phonon energy 
ui\ n in Eq. (24), listed in the table, exhibits an opposite 
trend to A. In fact, as represented in the MAD formula, 



the subtle balance of A and cu\ n determines T c 



F. Screened electron-electron interaction 

We next consider the RPA screened electronic inter- 
action W(r, r') [Eq. ( 16 )] , which is related to JC cl by 
Eq. (11). We calculated it with various settings of the 



origin and direction. We chose the four points as the 
origin r': The transition-metal, nitrogen, chlorine, and 
lithium sites. For the calculated directions, we selected 
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FIG. 9: (Color online) Our calculated RPA screened electronic interaction W(r,r') in Eq. (16 1 for /3-Lio.3HfNCl. The results 
are shown for the four different origins r'; r'=Rnf in the panel (a), t'=Rn in (b), r'=Rci in (c), and r'=RLi in (d), where 
Rx is the atomic coordinate of the species X. The calculations were performed along the two directions; one is the direction 
along the interlayer c axis, plotted as red crosses, and the other is taken to the parallel direction to an Hf-N bond in the ab 
plane as blue open circles. The geometrical configuration from the side view about the origin and the calculated direction is 
given in the inset. The figures also include the interaction based on TFA in Eq. (13 1 (solid lines). We do not show the values 
r— r'|<l (a.u.), because this regime is beyond the real-space resolution with the present plane-wave cutoff employed in the 



for 

screened-interaction calculation (12.8 Ry) 



80 



the following two: One is the direction along the inter- 
layer c axis and the other is taken to be parallel to an 
M-N bond in the ab plane. We show in Fig. [9] our cal- 
culated W(r, r') for /3-Lio.3HfNCl. The atomic geometry 
is depicted in the inset, where the origin and calculated 
direction are superposed as the green arrow. From the 
comparison among the four panels, we see significantly 
weak interactions around the Hf and N atoms building 
the block layer compared to the TFA one [solid lines, 
Eq. ( 13 )] , especially for the long-range part (>2 a.u.). In 



contrast, the RPA interaction around the CI and Li atoms 
seems to follow the TFA interaction. The strong screen- 
ing around the Hf and N atoms is because of the elec- 
tronic states near the conduction-band bottom formed 
by the Hf 5d and N 2p orbitalsP^ES T n the vicinity of r', 
the interaction around the Ni site seems to be smaller 
than that around the Hf one. This might be a conse- 



quence of the difference in the electronegativity (1.30 for 
Hf and 3.04 for N), but it should be noticed that this 
region is at the border of the resolution limit in the real 
space with the plane- wave cutoff of 12.8 Ry for screened 
Coulomb interaction in Eq. ( 17)p3The dependence of the 



RPA screened Coulomb interaction on the atomic species 
at the origin r' and the difference in the short-range part 
in the RPA and TFA interactions indicates that inhomo- 
geneity not described by TFA is essential for this system. 
To show this point clearer, we calculated W(r, r') for a 
simple metal Nb, shown in Fig. 10 The r' points were 



set to (0,0,0) (Nb site,'!" in the inset) and (1/2, 1/2, 
0) ("2"). The calculation was performed along the (111) 
direction (green arrow). We see that the W(r,r') hardly 
depend on the origin and the calculated directions, which 
is ascribed to a nearly uniform electronic density in the 
system. The interaction value itself is, contrary to the 
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cally follow those of //. However, as is well appreciated, 



the renormalization formula in Eq. ( 25 1 includes an em- 



r-r' [a.u.] 
FIG. 10: (Color online) Our calculated RPA screened elec- 



tronic interaction W(r, r') in Eq. ( 16 1 for Nb. The results are 



shown for the two different origins r'; one is set to the Nb 
atom (denoted by "1") and the other is on the center of the 
(100) surface ("2"). The calculations were performed along 
the (111) direction. The geometrical configuration is schemat- 
ically given in the inset, where purple spheres represent the 
Nb atoms and the direction taken in the W(r,r') calculation 
is indicated by green arrows. T he figures also include the 
interaction based on TFA in Eq. (13 1 (solid lines). 



/3-Li .3rIfNCl case, larger than the TFA one. 

The above-mentioned aspects of W(r,r') for j3- 
Lio.3HfNCl are common in essence to the other (3- 
Li^MNCl compounds. On the other hand, as the 
transition-metal element is replaced by a lighter one, we 
find that the reduction of the RPA screened interaction 
around the transition-metal site from the TFA one be- 
comes more appreciable. This may be due to the fol- 
lowing fact: When transition metal element is replaced 
from 5d to Ad or 3d, in general, the ionic radius becomes 
smaller and the electronegativity becomes larger. 

Next, to see the strength of the effective interaction 
between the electrons for min g each Cooper pair, we cal- 
culated n following Ea. Q, with /C cl - TFA in Eq. (l4|) 
and JC cl ' KPA in Eq. (20). The results for each material 
and each doping rate are given in Table |IV| As basic 
trends, we see that (i) /xrpa is considerably small com- 
pared to /iTFA (ii) the large doping rate gives the large /i 
value, and (hi) in the order of Hf, Zr, and Ti compounds 
(i.e., as the transition metal atom becomes lighter), the 
resulting fj, becomes larger. The trend (i) results from 
the damped RPA interaction around the nitrogen and 
transition-metal atoms. The trends (ii) and (iii) originate 
from the enhancement in the N(0) factor; although the 
averaged interaction ((/C c1 ))fs itself in Eq. (27) is weak- 
ened by the doping or almost unchanged with the sub- 
stitution from heavier to lighter element, the increase of 
N(0) (see Fig. [5] and Table IV) dominates the overall 
trend in \x over the ((/C o1 ))fs factor. We also give in the 
table the renormalized value \i* . The trends of /x* basi- 



pirical parameter E e \ and the setting detail of this param- 
eter strongly affect the final value. Hence, the estimate 
itself has little quantitative reliability. 



G. Transition temperatures 



We solved the SCDFT equation [Eq. Q or pB] with 
the constructed exchange-correlation kernels Z^, /C ph , 
and JC el from first principles. Figure [II] plots our calcu- 
lated superconducting gap as a function of temperature. 
The upper two panels (a) and (b) describe results for the 
Zr compounds and the lower two (c) and (d) describe 
those for the Hf compounds. The panels (a) and (c) [(b) 
and (d)] are the results where TFA (RPA) is adopted for 
/C cl . The pink (lighter) and blue (darker) squares in the 
figures represent the results for x=0.3 and 0.5, respec- 
tively, and the horizontal bars indicate the experimental 
T c range. 

We find two discernible differences between TFA and 
RPA, which are common in both the Zr and Hf com- 
pounds. First, unlike the results for Al and Nb in Sec. 
Ill, the estimated T c is higher in RPA than TFA. This is 
ascribed to the weaker RPA electronic interaction than 



the TFA one (Sec. |IVF). Second, while both TFA and 
RPA predict that an electron doping raises T c , the in- 
creasing trend of T c in RPA is less noticeable than that in 
TFA. This is because that the weakening of the RPA in- 
teraction in the MN layers (Fig.[9| compared to the TFA 
one becomes more appreciable in the low-doping regime 
than the high-doping one (as seen in Table IV ) , thus lead- 
ing to the relatively higher RPA T c than the TFA T c at 
x=0.3. The observed differences between TFA and RPA 
highlight a significance of the effects beyond TFA in (3- 
Li x MNCl. 

Importantly, our estimated T c is a few times lower 
than the experimental values, being in a sharp con- 
trast with the case of simple metals for which SCDFT 
successfully reproduces experimental T c . The disagree- 
ment observed for f3-Li x MNCl suggests that the present 
exchange-correlation kernels Z ph , /C ph , and /C cl is miss- 
ing something crucial to describe the superconductivity 
in this system. 

We summarize the SCDFT results in Table EJ We 
compare three T c s obtained from different treatments of 
the Coulomb interaction effect; ySCDFT-zc =o j or w } 1 j c } 1 
/C cl is neglected in the gap function, ySCDFT-TFA ca j_ 
culated with TFA, and T C SCDFT - RPA obtained by using 
RPA. For the Ti compounds, we found no superconduct- 
ing solution in calculations with JC cl . We also list T c es- 
timated by the MAD formula. The trend is basically the 
same as that of SCDFT. While the MAD gives higher 
T c , it should be noted that the MAD formula contains 
the empirical parameter. As a reference, we estimated 
A'* = Mscdft-rpa which reproduces our jjSCDFT-rpa f or 
the Zr and Hf compounds and list them at the bottom 
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FIG. 11: (Color online) Temperature dependence of the gap 
functions for /3-L^MNCl (cc=0.3 and 0.5; M=Zr and Hf). The 
panels (a) and (c) are the results with the Thomas-Fermi 
approximation for the electronic interaction kernel AC cl,TF in 
Eq. ( 14 1 , and (b) and (d) corres pond to the random-phase 



approximation for /C el,RP in Eq. (j20b . The experimental T c 



ranges are also shown. See the caption of Table [V] 



of the table. While the empirical formula [Eq. ( 25 )] gives 



TABLE V: Our calculated superconducting transition tem- 
peratures (K) with density functional theory for supercon- 
ductors and the McMillan- Allen-Dynes formula in Eq. ( 22 1 , 
and the parameter fi* given so that Eq. ( 22 1 reproduce 

/tiSCDFT-RPA ' — ' 



M 



Ti 



Zr 



Hf 



0.3 



0.5 



T SCDFT-/C-=0 ^ YJ 



0.3 0.5 0.3 0.5 



■* c 

thSCDFT-TFA 
■*- c 
thSCDFT-RPA 

Tf AU ( M *=0) 

71MAD / * \ 

J c IMtfaJ 

71MAD / * \ 



3.6 10.1 
<10" 3 <10~ s 



17.8 
0.8 
3.9 

18.8 
2.1 
6.2 



33.0 
2.2 
4.3 

31.2 
9.7 

11.5 



28.6 
4.4 
9.6 

27.1 
7.2 

14.8 



36.6 
7.2 
10.5 
26.9 
14.6 
16.6 



rr-iCXpt 
* 

MSCDFT-RPA 



10.0-15.2" 19.0-25.5" 
0.144 0.296 0.170 0.242 



0.113-0.196 (see Table IV), Mscdft-rpa is 0.144-0.296, 



" Ref. OH (/3-Li x ZrNCl) ; Ref . [21] (,3-Li x ZrNCl) ; 

Ref. [81] [^-Li x M„ ZrNCl; M= tetrahydrofuran (THF), 

propylene carbonate (PC)]; 
Ref. [82] (/3-A x ZrNCl; A=Li, Na, K); 
Ref. 1551 \8-LUM„ ZrNCl; M=(N, N)-dimethylformamide, 

dicmthylsulfoxidc] . 
b Ref. HI [0-Li o .48(THF)„HfNCl]; Ref. [72 (,3-Na x HfNCl); 
Ref.|2Q|(/9-Li x M H HfNCl; M=NH 3 , THF, PC); 
Ref.|8T](/3-A x M H HfNCl; A= Li, Na; M= THF, PC). 



which is much larger than the standard values (~0.10) 
used in literatures.^ These results suggest that suppres- 
sion of the Coulomb interaction by screening and renor- 
malization in the /3-MNC1 is not so strong. 



V. DISCUSSION 

In discussing overall aspect of the present SCDFT cal- 
culation, an analysis with the parameters in the MAD 
formula gives a simple insight because each parameter 
represents an average of the matrix quantity in SCDFT. 
We summarize in Fig.[l2]these parameters; A in the panel 
(a), Wln in (b), /i RPA in (c), and T C SCDFT " RPA in (d). The 
figure compares the three compounds (Ti, Zr, and Hf) 
and the empty and solid bars represent the results for 
x=0.3 and 0.5, respectively. While the material depen- 
dence of T c basically follows the A trend, the quantita- 
tive aspect of T c is determined as a balance in A, oj\ n , 
and \x. Because of the rather large /xrpa, the Ti com- 
pounds do not exhibit superconductivity. The doping 
dependence of T c also follows the A trend. However, be- 
cause of the decreasing (increasing) trend of U\ n (^rpa) 
upon the doping, the T c difference in x=0.3 and 0.5 is 
not substantial. Besides the overall trend of each factor, 
there is a notable quantitative disagreement in T c be- 
tween the theory and the experiment. Since the SCDFT 
assures a high accuracy of T c for a wide class of phonon- 
mediated superconductors, this disagreement is remark- 
able and should be seriously and carefully examined. To 
study this point in more detail, we calculated the doping 
dependence of T c fo r low er doping regime down to £=0.1. 
In the experiments, 19 20 T c is known to rapidly increase 
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as x is reduced down to ~0.05 for the Zr compounds, 
while it barely changes for £>0.1 for the Hf compounds. 
In order to examine whether this behavior is reproduced 
theoretically, we estimated T c for £=0.1 by SCDFT. In 
this calculation, the RPA electronic kernel in Eq. (20) 



is made with the rigid-band-shift approximation to the 
electronic state for £=0.3. For a 2 F, we employed that 
of £=0.3. This is justified from the fact that the DOS 
at Ep for £=0.1 is nearly the same as that for £=0.3 
(see Fig. [5]). The results are plotted in Fig. 13 and com- 
pared with the experiments. We see that the resulting 
T c lowers for small x, which is an opposite trend to the 
experiments. The doping dependence by SCDFT is rea- 
sonable because the screening due to the virtual states in 
the flat-band region just above Ep becomes less effective 
for small £ (as Ep is apart from this flat band). Indeed, 
the bare // is estimated to be 0.203 for the Zr compound 
and 0.237 for the Hf compound, which are respectively 



400 



a- 0.5 



2-10.0 




FIG. 12: (Color online) Transition-metal and doping depen- 
dence of the ab initio McMillan- Allen-Dynes parameters; A in 
the panel (a), cj ln in (b), ^ RPA in (c), and ySCDFT-RPA in ^ 
See also Tables |IV| and [V] Empty and filled boxes represent 
the values for the a;=0.3 and :r=0.5 compounds, respectively. 
For the panels (a)-(c), the boxes drawn in blue (red) indi- 
cate that the corresponding parameter contributes to raising 
(lowering) T c . 
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FIG. 13: (Color online) Doping dependence of transition tem- 
peratures derived from SCDFT in /J-Li^ZrNCl (upper panel) 
and /3-LbTlfNCl (lower panel). The errors due to the random 
sampling in the SCDFT calculation are also shown with the 
vertical bars. Experimental transition temperatures are taken 
from Ref ll9l for the Zr compound and Ref l2UI for the Hf one. 
The notation THF stands for tetrahydrofuran. The dashed 
lines are guides to the eyes. 



larger than those for £=0.3 (0.161 for both of the Zr and 
Hf compounds). 

The present SCDFT is based on the standard ME the- 
ory, in which one neglects the momentum dependence of 
the quantities. Also, the screened Coulomb interaction 
is described with RPA and the frequency dependence is 
ignored. The present discernible difference in T c or in 
its doping dependence between the theory and experi- 
ment suggest the importance of the effect not included 
in the analysis in the ME level. In fact, in the optical 
experiment,^ /3-LizZrNCl is known to have the plasma 
edge in the low-energy Drude band less than ~0.1 eV, 
exhibiting a redshift with the lowering doping amount. 
Since the presence of the low-energy plasma edge af- 
fects the screened Coulomb interaction in low-frequency 
regime, it can be relevant to the low-energy properties of 
the present system as proposed earlier in Ref. 40, Mean- 
while, there remains a possibility that correction in the 
Kohn-Sham orbitals improves the present error; for ex- 
ample, the use of the hybrid functional in deformation- 
potential calculation seems to lead to the increase of the 



15 



electron-phonon couplingPSlt will be interesting that the 
present SCDFT calculation is performed with such hy- 
brid functional to study the degree of the improvement. 



calculations reasonably reproduces the experimental 
results, as well as the previous theoretical results. 



VI. SUMMARY AND OUTLOOK 



We have performed a comprehensive SCDFT analysis 
for doped /3-MNC1, where M=Ti, Zr and Hf, to study the 
effects described within the Migdal-Eliashberg theory 
without introducing any adjustable parameters. On the 
basis of the ab initio electronic band structure, phonon 
spectrum, electron-phonon coupling, and electron- 
electron interaction in the RPA level, we estimated 
superconducting T c . The resulting T c values are <4 
K for the Zr compounds and <10 K for the Hf ones, 
whereas the Ti compounds do not show superconduc- 
tivity because of the large electronic interaction. While 
the dependence of T c on the transition-metal species 
is consistent with the experiments, our calculated T c 
values are less than a half of the experimental values 
(10-15 K for the Zr case and 19-26 K for the Hf one) 
and the theoretical doping dependence of T c is opposite 
to that of the experiment. Although we do not rule 
out some possibilities of improvement by corrections to 
one-particle states as the starting point, the revealed 
discrepancy between the theory and experiment, as is 
known in the case of the unconventional superconduc- 
tors, appeals the relevance of something ignored in the 
standard Migdal-Eliashberg theory to the superconduct- 
ing mechanism of the doped /3-MNC1. This remains an 
open question. 



TABLE VI: Our calculated frequencies (cm 1 ) of the Raman- 
active modes for undoped /3-MNC1. 



M 


Ti 




Zr 






Hf 






present 


present 


cxpt. a 


theory 


present 


cxpt. a 


thcory c 


Ai s #l 


603 


598 


591 


564 


634 


612 


604 


Ai s #2 


391 


333 


331 


312 


205 


292 


386 


Ai s #3 


214 


196 


191 


182 


169 


163 


191 


E,#l 


640 


622 


605 


020 


652 


033 


- 


E,#2 


227 


191 


184 


171 


166 


156 


- 


E,#3 


133 


130 


128 


114 


113 


110 


- 



Ab initio DFT-LDA pseudo-potential calculation in Ref. I29l 
12 Ab initio DFT-LDA pscudopotcntial calculation in Ref. I32l 
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Appendix A: T-point Phonons 

To check the reliability for phonon properties, we 
calculated vibrational frequencies at the Y point and 
compared with previous theoretical and experimental 
results. Tables IVII and IVIII list the Raman and infrared 
frequencies for the undoped system, respectively. Also, 
Tables |VHI| and |IX| are the results for doped samples. 
For both of the undoped and doped samples, our 



TABLE VII: Our calculated frequencies (cm -1 ) for infrared- 
active modes for undoped /3-MNC1. The values in the paren- 
theses are the frequencies of the longitudinal vibrations. 



M 


Ti 




Zr 




Hf 






present 


present 


cxpt. a 


theory 


present 


thcory c 


A 2 „#l 
A 2u #2 
E„#l 
E„#2 


638(666) 
303(317) 
537(736) 
184(206) 


630(673) 
268(286) 
538(713) 
171(198) 


529 

666? 

165 


597(630) 
257(272) 
547(680) 
155(179) 


638(684) 
257(276) 
565(721) 
159(189) 


617 
128 
632 
105 



a Infrared spectroscopy taken from Ref. |72| 

Ab initio DFT-LDA pscudopotcntial calculation in Ref. |29l 
c Force-constant-modcl calculation reproducing 
the Raman spectra in Ref. 1851 
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TABLE VIII: Our calculated frequencies (cm ) of the 
Raman-active modes for doped /3-Li^MNCl. 



M 


Ti 




Zr 




HI 






present 


present 


expt. a 


theory" 


present 


expt. c 


X 


0.3 0.5 


0.3 0.5 


0.16 


1/6 


0.3 0.5 


0.35 


Ai 8 #l 


594 556 


557 473 


582 


565-548 


589 496 


595 


Ai s #2 


369 372 


304 291 


322 


303-296 


250 228 


282 


Ai s #3 


210 223 


198 209 


188 


186-217 


181 197 


147 


E,#l 


634 633 


603 616 


608 


599-633 


607 631 


616 


E s #2 


222 224 


181 170 


178 


169-165 


138 134 


157 


E,#3 


125 116 


123 111 


123 


116-131 


119 110 


106 



Raman spectroscopy taken from Ref. |72| 
Ab initio DFT-LDA pseudopotcntial calculations in Ref. 1 29 1 
"Raman spectroscopy for /S-Na^HfNCl taken from Ref. 1851 



TABLE IX: Our calculated frequencies (cm 1 ) for infrared- 
active modes for doped /3-Li x MNCl. 



M 


Ti 


Zr 


Hf 






present 


present theory a 


present theory" 


X 


0.3 0.5 


0.3 0.5 1/6 


0.3 0.5 


0.5 


A 2u #l 


638 622 


625 585 612-583 


632 606 


596 


A 2 „#2 


283 346 


247 291 314 


254 300 


108 


A 2 „#3 


223 152 


198 142 235-184 


186 128 


183 


E u #l 


588 587 


557 540 549-603 


576 534 


618 


E u #2 


231 271 


215 221 262 


217 232 


85 


E u #3 


127 86 


119 83 148-113 


108 77 


146 



1 Ab initio DFT-LDA pseudopotcntial calculations 
in Ref. (29) 

Force-constant-modcl calculation reproducing 
the Raman spectra for /3-Nao.5sHfNCl in Ref. |85| 
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